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(57) Abstract: The invention relates to a 
method for the computer-aided reconstruction 
of a three-dimensional anatomical object (3) 
from diagnostic image data. First of all, a 
diagnostic image data set of the object (3) is 
acquired. Then a seed point (5) is set, starting 
from which the object is reconstructed within 
a reconstruction volume (4). Thereafter, an 
adjacent point of the reconstruction volume 
(4) likewise belonging to the object (3) is 
located in accordance with a propagation 
criterion, which is calculated by means 
of a mathematical analysis of local areas 
(6, 7), assigned to the point concerned, of 
the image data set. Reconstruction of the 
three-dimensional structure of the object (3) 
is then performed within the reconstruction 
volume (4) by multiple repetition of this 
method step and propagation along the 
located adjacent points. To apply such a 
reconstruction method to image data obtained 
by means of rotational X-ray imaging, wherein 
a plurality of two-dimensional projection 
images (1, 2) are recorded from different 
projection directions, the invention proposes 
that the propagation criterion be calculated 
by subjecting the local image areas (6, 7) of 
the two-dimensional projection images (1, 2) 
in each case individually to the mathematical 
analysis. 
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Method for the reconstruction of three-dimensional objects 



The invention relates to a method for the computer-aided reconstruction of a 
three-dimensional anatomical object from diagnostic image data, having the method steps: 

a) acquisition of a diagnostic image data set of an object, 

b) setting of a seed point belonging to the object within a reconstruction 

5 volume, 

c) location of an adjacent point, likewise belonging to the object, within the 
reconstruction volume in accordance with a propagation criterion, which is calculated by 
means of a mathematical analysis of local areas, assigned to the point concerned, of the 
image data set, 

10 d) reconstruction of the three-dimensional structure of the object within the 

reconstruction volume by multiple repetition of method step c) and propagation along the 
adj acent points thus located. 

In addition, the invention relates to a computer program and an imaging 
apparatus with* computer means for performing this method. 

15 In the field of angiography, three-dimensional medical imaging methods, such 

as for example three-dimensional rotational X-ray imaging (3D-RX) or magnetic resonance 
imaging (MRI), are growing in importance. The volume image data obtained with such 
methods contain interesting information for diagnosis of vessel diseases, such as for example 
stenoses or aneurysms. In such cases, visualization of the vessel structures is crucial in 

20 allowing a doctor treating the condition to recognize quickly and reliably potential danger 
sources (e.g. an impending infarction or thrombosis). 

Computer-aided three-dimensional reconstruction of the vessel system of a 
patient from the image data acquired on the one hand allows the profile of the blood vessels 
to be visualized with high reproduction accuracy, anatomical structures not belonging to the 

25 vessel system concerned being hidden. On the other hand, the three-dimensional 

reconstruction of the vessel structures is a useful aid in planning interventions, such as for 
example left coronary catheter investigations (PTC A). 
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A three-dimensional reconstruction method for analyzing volume image data 
acquired by magnetic resonance angiography (MRA) is known for example from an article 
by Young et al (S. Young, V. Pekar and J. Weese, „Vessel Segmentation for Visualization of 
MRA with Blood Pool Contrast Agent", MICCAI 2001, 491-498, Utrecht, October 2001). 
5 The previously known method serves, inter alia, to separate the arterial and venous vessel 

systems from one another during visualization of the image data. According to the previously 
known method, first of all a diagnostic image data set is acquired in the form of a volume 
image of the vessel structures of interest, using a suitable contrast agent. Then, a user sets a 
seed point within a reconstruction volume, this seed point being identified by the user as 

1 0 belonging to a venous vessel. Automatic three-dimensional reconstruction of the selected 

vessel then takes place by means of a propagation method, which is based on a mathematical 
analysis of the respective local image areas. Starting from the seed point, points within the 
reconstruction volume are identified, in accordance with a propagation criterion supplied by 
the mathematical analysis, as belonging or not belonging to the vessel, whereby segmentation 

15 of the reconstruction volume takes place. Propagation continues until the entire structure has 
been reconstructed or until a set end point is reached. The mathematical analysis applied for 
calculation of the propagation criterion is of fundamental importance to the previously known 
method. In the stated article, a mathematical filter is proposed in this respect, which is based 
on evaluation of the second derivatives of the gray scale values within the local image areas. 

20 A proposed alternative involves adaptation of the local image data to a cylinder model, by 
means of which the mathematical analysis is rendered selective for image structures typical 
of blood vessels. 

In rotational X-ray imaging, a plurality of two-dimensional projection images 
is recorded at different projection angles, for example by means of a C-arm X-ray apparatus. 

25 To make the blood vessels of the patient under investigation visible in the projection images, 
an X-ray absorbent contrast agent is injected into the patient. A problem with this 
investigation method is that the blood vessels typically have a complicated three-dimensional 
profile, which it is difficult for the doctor to detect solely on the basis of two-dimensional 
projection images. The missing three-dimensional information within a projection image 

30 must be added by the doctor by comparison with images recorded at other projection angles. 

It is now possible to generate a volume image data set from the plurality of 
two-dimensional projection images recorded by means of 3D-RX using suitable modeling or 
back projection methods on a suitable computer. This volume image data set may then 
undergo an analysis of the type outlined above for the purpose of reconstruction of the three- 
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dimensional vessel structures. This procedure, however, is disadvantageously associated with 
considerable computing power. A further disadvantage is that, in particular if the coronary 
vessels of the patient are to be investigated, generation of the projection images has to be 
ECG-controlled, so that the coronary arteries are recorded in all the images in the same phase 
5 of the heart beat cycle. Because of the need for ECG control, only a comparatively small 
number of images is then available for each phase of the heart beat cycle, which means that 
the volume images reconstructed therefrom reproduce the vessel structures only relatively 
inaccurately. A quantitative analysis according to the above-described reconstruction method 
does not then provide any usable results. 

10 

Taking this as basis, it is an object of the present invention to provide a 
method of segmenting a reconstruction volume which is in a position, starting from a 
comparatively small number of two-dimensional projection images, to determine the three- 
15 dimensional structure of the object reliably, precisely and using as little computing power as 
possible. 

In the case of a method of the above-mentioned type, this object is achieved 
according to the invention in that, in method step a), a plurality of two-dimensional 
projection images is recorded from different projection directions, the propagation criterion 

20 being calculated in method step c) by subjecting the local image areas of the two-dimensional 
projection images in each case individually to mathematical analysis. 

The basic concept of the invention is to perform the computer-aided 
segmentation of the reconstruction volume directly by means of a propagation method known 
per se, without any intermediate reconstruction of a three-dimensional volume image data set 

25 from the projection images. In the process, propagation in the reconstruction volume along 
the contours of the object to be reconstructed is controlled by combining the information 
obtained by means of the mathematical analysis applied to the individual two-dimensional 
projection images to yield a uniform propagation criterion. 

To this end, it is possible, for example, to identify a point in method step c) as 

30 belonging to the object, provided that the mathematical analysis yields a result which agrees 
for a plurality of two-dimensional projection images. This procedure takes account of the fact 
that, on the basis of projection, the mathematical analysis of an individual, two-dimensional 
projection image may cause the point concerned to appear to belong to the object even when 
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this is not actually the case. Only a comparison with the results obtained by mathematical 
analysis of the other projection images in relation to this point allows reliable segmentation. 

The local image areas are appropriately determined in method step c) by 
projecting the point concerned within the reconstruction volume in accordance with the 
5 respective projection directions into the image planes of the two-dimensional projection 
images. In this way, the geometric conditions when the projection images are recorded are 
replicated, in order to be able to achieve assignment of the points of the reconstruction 
volume and the image points of the two-dimensional projection images. 

By the mathematical analysis in method step c), a propagation coefficient 

10 ought appropriately to be calculated in each case as propagation criterion for each two- 
dimensional projection image, the value of which coefficient indicates whether the point 
concerned belongs to the object or not. Such a coefficient is particularly well suited to 
performance of the method according to the invention by means of a computer, since location 
of points belonging to the object to be reconstructed may be effected by simple numerical 

1 5 comparison. For example, the procedure may be performed in such a way that, in method 
step c), a point is identified as belonging to the object, provided that the propagation 
coefficient assumes a large value for a plurality of two-dimensional projection images. 

A characteristic of blood vessels is their axial symmetry. They extend a long 
way in one direction and only a short way in the direction perpendicular thereto. This 

20 morphological characteristic may be used according to the invention to calculate the 
propagation coefficient. For three-dimensional reconstruction of vessel structures, it is 
accordingly sensible, during calculation of the propagation coefficient, to calculate the 
inherent values of the Hesse matrix of the gray scale values in the local image area of the 
respective two-dimensional projection image. By evaluating these inherent values, 

25 propagation then follows the image structures with - from a spatial point of view - the lowest 
possible gray scale curvature values, because the Hesse matrix provides information about 
the local second derivatives of the gray scale values. Suitable formulae for calculating the 
propagation coefficient on the basis of the inherent values of the Hesse matrix may be found, 
for example, in the above-cited article by Young et al. The propagation coefficient may be 

30 calculated from the two-dimensional projection images for example as follows: 

0 A 2 (*)>0, 

A 2 (3c)<0, >« =pq> 5 = VEA 2 • 
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In this equation, a and c are weighting factors and X\ ( x ) and A,2 ( x ) are the 
inherent values of the local gray scale value Hesse matrix calculated at the point x within the 
respective two-dimensional projection image. More details about this may be found in the 
above-cited publication by Young et al. 
5 When calculating the propagation coefficient for the respective two- 

dimensional projection image, adaptation to a cylinder model within the local image area 
may also be calculated. Such a cylinder model, which is also described in detail in the stated 
article by Young et al, likewise makes vessel structures distinguishable from other 
anatomical structures. 

10 Reconstruction is appropriately stopped when a predeterminable end point is 

reached during propagation in method step d). Such an end point may either be 
predetermined interactively or determined automatically, for example on the basis of the size 
of the reconstruction volume. 

An imaging apparatus, in particular a C-arm X-ray apparatus, for performing 

15 the method according to the invention constitutes the subject matter of claim 9, according to 
which a computer means of the imaging apparatus is provided with a program such that the 
two-dimensional projection images are recorded according to the above-described method. 
For angiographic investigations of the coronary arteries, the imaging apparatus appropriately 
comprises ECG control as claimed in claim 10, so as to be able to record the projection 

20 images synchronously with the heart beat. 

A computer program as claimed in claim 1 1 is suitable for performing the 
method according to the invention, for example on an imaging apparatus equipped with a 
suitable computer means. The software required therefore may be made available to the users 
of corresponding imaging apparatus advantageously on a suitable data medium, such as a 

25 floppy disk or a CD-ROM, or by downloading from a data network (Internet). 



The invention will be further described with reference to examples of 
embodiments shown in the drawings to which, however, the invention is not restricted. In the 
30 Figures: 

Fig. 1 is a schematic representation of the method according to the invention 
for reconstructing a three-dimensional anatomical object; 

Fig. 2 shows an imaging apparatus according to the invention. 
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Fig. 1 shows a diagnostic image data set consisting of two two-dimensional 
projection images 1, 2, which image data set was acquired by means of X-ray fluoroscopy. 
Each of the projection images 1, 2, recorded at different projection angles, shows a branched 
5 blood vessel 3 of a patient. The projection images 1 5 2 accordingly show the same blood 
vessel 3 from different perspectives. To acquire the image data set, a contrast agent was 
administered to the patient, such that the blood vessel 3 shows up dark in the projection 
images. To reconstruct the three-dimensional structure of the blood vessel 3 according to the 
invention, a seed point 5 is firstly set within a reconstruction volume 4. The contour of the 

10 blood vessel 3 is then reconstructed in the volume 4, by locating adjacent points in the 
volume 4 in each case belonging to the blood vessel 3 in accordance with a propagation 
criterion. To this end, local image areas 6 and 7 belonging to the respective point 5 within the 
two-dimensional projection images 1 and 2 respectively are in each case subjected 
individually to mathematical analysis. After location of a point adjacent to the seed point 5, 

15 the procedure is repeated for points in turn adjacent to this point, until the entire structure of 
the blood vessel 3 has been reconstructed within the volume 4. The point investigated in each 
case with each propagation step" is fdehtified~as belonging to the blood vessel if the 
mathematical analysis of the local image areas 6 and 7 gives a positive result for both 
projection images 1 and 2 respectively. The local image areas 6 and 7 are determined by 

20 projecting the point 5, in accordance with the projection directions in which the two images 1 
and 2 were recorded, into the image planes of these two images. This is indicated in Fig. 1 by 
arrows 8 and 9. 

The imaging apparatus illustrated in Fig. 2 is a C-arm X-ray apparatus, which 
comprises a C-arm 10, which is suspended by means of a holder 1 1 from a ceiling (not 

25 described in any more detail). An X-ray source 12 and an X-ray image converter 13 are 

guided movably on the C-arm 10, such that a plurality of two-dimensional projection X-ray 
images of a patient 15 lying on a table 14 in the center of the C-arm 10 may be recorded at 
different projection angles. Synchronous movement of the X-ray source 12 and the X-ray 
image converter 13 is controlled by a control unit 16. During image recording, the X-ray 

30 source 12 and the X-ray image converter 13 travel synchronously around the patient 15. The 
image signals generated by the X-ray image converter 13 are transmitted to a controlled 
image processing unit 17. The heart beat of the patient 15 is monitored using an ECG 
apparatus 18. The ECG apparatus 18 transmits control signals to the image processing unit 
17, such that the latter is in a position to store a plurality of two-dimensional projection 
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images in each case in the same phase of the heart beat cycle, in order in this manner to 
perform an angiographic investigation of the coronary arteries. The image processing unit 17 
comprises a program control, by means of which three-dimensional reconstruction of a blood 
vessel detected with the image data set thus acquired is performed, according to the above- 
5 described method. The reconstructed blood vessel may then be visualized in known manner 
on a monitor 19 connected to the image processing unit 17. 
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CLAIMS: 



1 . A method for the computer-aided reconstruction of a three-dimensional 
anatomical object (3) from diagnostic image data, having the method steps: 

a) acquisition of a diagnostic image data set of the object (3), 

b) setting of a seed point (5) belonging to the object (3) within a reconstruction 

5 volume (4), 

c) location of an adjacent point, likewise belonging to the object (3), within the 
reconstruction volume (4) in accordance with a propagation criterion, which is calculated by 
means of a mathematical analysis of local image areas (6 5 7), assigned to the point (5) 
concerned, of the image data set, 

10 d) reconstruction of the three-dimensional structure of the object (3) within the 

reconstruction volume (4) by multiple repetition of method step c) and propagation along the 

adjacent points thus located, 

characterized in that, in method step a), a plurality of two-dimensional 

projection images (1, 2) is recorded from different projection directions, the propagation 
15 criterion being calculated in method step c) by subjecting the local image areas (6, 7) of the 

two-dimensional projection images (1, 2) in each case individually to the mathematical 

analysis. 

2. A method as claimed in claim 1, characterized in that, in method step c), a 
20 point is identified as belonging to the object (3) if the mathematical analysis yields a result 

which agrees for a plurality of the two-dimensional projection images (1, 2). 

3. A method as claimed in one of claims 1 or 2, characterized in that, in method 
step c), the local image areas (6, 7) are determined by projecting the respective point (5) 

25 within the reconstruction volume (4) in accordance with the respective projection directions 
into the image planes of the two-dimensional projection images (1,2). 



4. A method as claimed in one of claims 1 to 3, characterized in that, a 

propagation coefficient is calculated by the mathematical analysis in method step c), as a 
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propagation criterion for each two-dimensional projection image (1, 2), the value of which 
coefficient indicates whether the point (5) concerned belongs to the object or not. 

5. A method as claimed in claim 4, characterized in that, during calculation of the 
5 propagation coefficient, the inherent values are calculated of the Hesse matrix of the gray 

scale values in the local image area (6, 7) of the respective two-dimensional projection image 
(1,2) 

6. A method as claimed in claim 4, characterized in that, when calculating the 
10 propagation coefficient for the respective two-dimensional projection image (1, 2), an 

adaptation to a cylinder model within the local image area (6, 7) is calculated. 

7. A method as claimed in one of claims 4 to 6, characterized in that a point is 
identified in method step c) as belonging to the object (3) if the propagation coefficient 

15 assumes a large value for a plurality of two-dimensional projection images (1,2). 

8. A method as claimed in one of claims 1 to 7, characterized in that the 
reconstruction is stopped when a predeterminable end point is reached during propagation in 
method step d). 

20 

9. An imaging apparatus, in particular a C-arm X-ray apparatus, having means 
(10, 11, 12, 13, 16) for generating an image data set, which set comprises a plurality of two- 
dimensional projection images of a body part of a patient (15) recorded from different 
projection directions, and having computer means (17) for reconstructing a three-dimensional 

25 anatomical object from the image data set, 

characterized in that the computer means (17) comprise a program control 
which operates according to the method as claimed in one of claims 1 to 8 to reconstruct the 
object. 

30 10. An imaging apparatus as claimed in claim 9, characterized by an ECG control 

(18), by means of which recording of the two-dimensional projection images can be 
controlled in accordance with the heart beat cycle of the patient (15). 
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11. A computer program for an imaging apparatus as claimed in claim 9 or claim 

10 5 characterized in that the method as claimed in one of claims 1 to 8 is implemented by the 
computer program on the computer means of the imaging apparatus. 
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